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Abstract 

This article concerns optimal estimates for non- homogeneous degenerate elliptic equa- 
tion with source functions in borderline spaces of integrability. We deliver sharp Holder 
continuity estimates for solutions to p-dcgcncratc elliptic equations in rough media with 

sources in the weak Lebesgue space i^^ak- ^^e borderline case, / € -^Xcak' solutions 
may not be bounded; nevertheless we show that solutions have bounded mean oscillation, 
in particular John-Nirenberg's exponential integrability estimates can be employed. All 
the results presented in this paper are optimal. Our approach is inspired by a power- 
ful Caffarelli-type compactness method and it can be employed in a number of other 
situations. 
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1 Introduction 

Central theme in the theory of elliptic partial differential equations, the classical Poisson 
equation 

-Au = f{X), (1.1) 

models important problems from theoretical physics, mechanical engineering to biology, eco- 
nomics, among many other applications. One of the key objectives in the analysis of Poisson 
equations is to assure regularity of u based on smoothness or integrability properties of its 



laplacian, /. In this context, Schauder estimates is a fundamental result. It assures that the 
Hessian of u, D^u, is as regular as /, provided / has an appropriate modulus of continuity. 

h/2) 



More precisely, if / G C"(Si), < a < 1 then u G C^'''{Bi/2) and 



lkllc2>-(Bi/2) ^ C"" {ll/llc-(Bi) + j (1-2) 

for a dimensional constant C„. Schauder estimate is sharp in several ways. Clearly if n G 
C^'", then its laplacian is a-Holder continuous. Also if / is merely continuous, one cannot 
assure u G C^, nor even C^^l bounds are available. Schauder estimates also fail in the upper 
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extreme, a = 1, i.e., if / G Lip, it is not true in general that u G Cj^'^. 

Establishing regularity of solutions to (jl.ip reduces to understanding the behavior of the 
Newtonian potential of /, 

^/W ■■= J Jx^Yf^f^^^dY. (1.3) 

The kernel that appears in ()1.3|) . r(X) = \X — is the fundamental solution of the 

laplacian. The second derivative of F, DijT ~ |X — Y\~"' is not integrable, but it is almost 
integrable, in the sense that \X — Y\'^DijT is integrable for any < e. This is the key 
observation that explains why Schauder estimates hold when / G C", < a < 1, and it fails 
when / is merely bounded or continuous. 

In several applications, the source function / is not continuous, but only (/-integrable, 
i.e., / G L''{Bi), for some 1 < g < oo. In this case, the corresponding regularity theory, due 
to Calderon and Zygmund, asserts that u G W^''^{Bi^2) 

lkllvi/2>9(Bi/2) < Cn{\\f\\Li[Bi) + WuhiiBi)} , (1-4) 

In particular, if / G then u G W'^''^ for all g < oo and by Sobolev embedding, 
u G C^'°' for any a < 1. This type of thesis is usually called almost optimal regularity 
result. Heuristically, for borderline hypotheses, almost optimal regularity result is the best 
one should hope for. 

Regularity theory for problems in rough heterogeneous media, i.e., when governed by 
elliptic equations with measurable coefficients, is rather more sophisticated, and even for the 
homogeneous equation 

V • {aij{X)Du) = 0, 

solutions are, in general, known to be only Holder continuous. This is the content of De 
Giorgi, Moser and Nash regularity theory. C alder on- Zygmung regularity estimates are not 
available in this setting. In even more complex models, the laplacian in (jl.ip is replaced by 
further involved nonlinear elliptic operators, 

-V ■a{X,Du) = f{X), (1.5) 

where a: Bi x R" — >■ M" is p-degenerate elliptic vector field. Throughout this paper we shall 
always assume the following standard structural assumption on the vector field a: 

f \a{x,o\ + \d^a{x,om < Aieri 
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for positive constants < A < A < +00. As usual in the literature, we could also include 
a parameter s > as to distinguish the model p-Laplacian operator (s = 0) from the 
nondegenerate one (s > 0), see for instance [14:\ I15|. Throughout this paper, constants that 
depend only upon n, p, A and A will be called universal. 

We recall that Equation (II. 5p appears for instance as the Euler-Lagrange equation of the 
minimization problem 

j F{X, Vu) + f{X)udX min, 

where the variational kernel F{X, is convex in ^, F{X, ~ 1^1^ and F{X, X^) = \X\pF{X, 
A typical operator to keep in mind is the p-laplace in rough media, 

-V ■ {aij{X)\Vu\P~^Vu), 

where bounded, positive definite matrix. 

The regularity theory for Equation (jl.Sp is nowadays fairly well established; however 
it is considerably more subtle than the corresponding linear, uniform elliptic theory. For 
instance, it is well known that p-harmonic functions are locally C^'" for some a that depends 
on dimension and p. The precise value of optimal a is, in general, unknown. 

The main goal of the present article is to determine optimal and almost optimal regularity 
estimates for solutions to Equation (jl.Sp , based on integrability properties (or more generally 
on the behavior of the distributional function) of the source /. The regularity estimates 
presented in this paper do not depend much on the concept of weak solution used. Indeed, 
they can be understood as a priori estimates that do not depend on any further regularity 
property of / or u. In the proofs, though, we shall always work with distributional solutions. 
However the same arguments go through, with no change, if one chooses to use the notion of 
entropy solutions, see [4] or any appropriate approximation scheme. Also we mention that, 
per our primary motivations, throughout the whole paper we shall only consider the range 

l<p<n, (1.7) 

where n is the dimension and p is the degeneracy exponent of the vector field a{X, Du). 

For L°°-bounds of solutions to Equation (jl.Sp . the borderline integrability condition on 
the source function / is . More precisely, if / G Lp , for any tin^ e > 0, solutions 
are bounded; however one cannot bound the L°°-norm of u by the Lp norm of /. The 
first result we show in this paper. Theorem 13. H is an optimal BMO estimate of solutions 

with source functions in the weak Lebesgue space L^^^^. Under slightly different structure 
assumptions, a similar result has been obtained by G. Mingione, Theorem 1.12 in [T3], as a 
consequence of potential analysis considerations (see also [iHj). Our proof is neither based 
on potential analysis nor on singular integral considerations. Instead, it is inspired by a 
powerful compactness type of argument, see [5], [6], and also [U [2]. The case p = n, i.e., 
for the n-Laplacian equation, with / being a finite measure relates to the article [10]. These 
results could be delivered by our methods as well. We emphasize that in the case p = n, 
^weak functions may not define a finite measure. Nevertheless, Theorem 13.11 provides a priori 
estimates for a priori regular solutions. When / is also a measure then this implies an 
existence and regularity theorem together with known approximation machineries. 

As soon as the source function / becomes (^ + e)-integrable, we show that solutions are 
in fact continuous. Not only do we show continuity of solution, but actually we provide the 
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precise sharp Holder exponent of continuity of u based only on the integrability of / and the 
regularity theory available for a-harmonic functions. Once more, the proof of such a result is 
based on compactness method and explores only the behavior of the distributional function 

of the source /, that is, / needs only to belong to the weak Lebesgue space i^J^^k' 1 < ^ < P- 
In this case, we show, see Theorem 14.11 and Theorem 14.21 that 

min^ — H— -an t 

^ ^ ^loc 

where ao is the universal optimal Holder exponent for solutions to —V • a{X, Du) = 0. 
Furthermore, we obtain the appropriate a priori universal estimate. Such a result brings 
important novelties. The first one is the optimal regularity space u lies. In many applications, 
for instance in free boundary and geometric problems, it is important to determine accurately 
how fast the solution grow away from its zero level set. In such a setting, knowing the precise 
regularity estimate is crucial for the program. Example of such problems are equations with 
singular terms, —ApU ~ u~'^ , 7 > 1. For these free boundary geometric problems, solutions 
are expected to behave like l^l'', near a free boundary point. Thus it is important to establish 

regularity estimates where potentials are assumed to belong to L^'^^^, but not in the classical 

n 

Lebesgue space L/'t. Another important advantage of our approach concerns its flexibility, 

which allows further generalizations, for instance to equations with measure data, to systems, 

or even to p-degenerate equations in nondivergence form, J-{X, u, Vu, D'^u) = f, where 

^s, ^p, ^M) ~ (,P~^T{X, s,p, M), for ,^ > 0. For this class of problems, compactness is 

consequence of Harnack type inequality as in the original approach in [5]. When projected 

to the constant coefficient case, the optimal estimate established in this paper is in 

accordance to the gradient estimates obtained in [15^19] through a powerful and sophisticated 

011 

nonlinear potential theory. Indeed, for the model equation —ApU = / G LJ^^^, it follows 

9n(p-l) 

from [15], Theorem 1 and p]. Theorem 1.1 that Vn G -^wcak* ' ^^'^^ by Morrey embedding 

p 9-1 

Theorem, u G Cp-^' * . 

The paper is organized as follows. In Section [2] we prove a basic compactness Lemma 
which assures that if ||/|| ^ is small, then there exists an ao-Holder continuous function 

"^wcak 

close to u in LP{Bi/2)- Section [3] is devoted to the proof BMO estimates. In Section U] we 
address the optimal C" regularity theory. 

Acknowledgement. The author would like to express his gratitude to Giuseppe Mingione, 
for several insightful comments that benefited a lot the final presentation of this article. The 
author also thanks the anonymous referee for such a careful revision. This work is partially 
supported by CNPq-Brazil. 



2 Compactness of solutions 

In this section, we establish a compactness result for solutions to non-homogeneous Poisson 
Equations (jl.Sp that will play a fundamental role in the proof of Theorem 13.11 and Theorem 
14. 1[ In fact Lemma l2.ll follows as a consequence of Lemma 3.2 in [8]. We include here a 
proof for completeness purposes. 



4 



Lemma 2.1. Let u G W^'^{Bi) be a weak solution to (|1.5p . with f \uY'dX < 1. Given 

6 > 0, there exists a < e ^ 1, depending only on p, n, X, A, u and 6, such that if 

WfWj .^<s, (2.1) 

then there exists a function h in B1/2 satisfying 

- V • a(X, V/i) = 0, inBi/2, (2.2) 
for some vector field a satisfying (jl.6p with the same ellipticity constants A and A such that 

\u{X)-h{X)\PdX <5P. 



Proof. Let us assume, for the purpose of contradiction, that the thesis of the Lemma fails. 
If so, there would exist a (5o > and sequences 

«fe G VFi'P(i?i), and /fcG4eak(^i)' 

satisfying 

/ \uk{X)Y'dX <l, (2.3) 

JBi 

for all /c > 1, 

-V-afc(X,Vnfc) = A inSi, (2.4) 
where Ok satisfies ()1.6p with constants A and A and 

ll/fcll ^ =0(1), (2.5) 

weak ^ ^ ' 

as A; —7- 0; however 

|n,.(X)-/i(X)|PdX>,5^, (2.6) 



/ 



for any solution h to the homogeneous problem ()2.2p in i?i/2 and all /c > 1. 

Now by standard Caccioppoli's type energy estimates, see for instance. Theorem 6.5 and 
Theorem 6.1 in |llj (notice that < ^ within the range 1 < p < n), we verify that there 
exists a constant C = C{n, A, A) such that 

\Vuk\^dX < C, 

Bl/2 

for all k > 1. Thus, up to a subsequence, there exists a function u G W^'P{Bi/2) for which 
Uk u mW^'P{Bi/2) and Uk ^ u in LP {B1/2). (2.7) 

In addition, in view of (|2.4p and (|2.5p . by classical truncation arguments, see for instance [3], 
we know 

Vnfc(X) -)■ Vu{X) for a.e. X G B^/2. (2.8) 
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Furthermore, by Ascoli Theorem, up to a subsequence, the sequence of vector fields afc(X, •) 
converges locally uniformly to a vector field a satisfying ()1.6p . Given a test function (p £ 
WQ'^iBi/2), from ([23]), ([22]) and (USD we have 

/ a{X,Du) -VcpdX = I ak{X,Duk)-V<l)dX + o{l) 

= / h^dX + o{l) 
= o(l), 

as fc — )• 00. Since (f) was arbitrary, we conclude u is a solution to the homogeneous equation 
in i?i/2- Finally we reach a contradiction in (j2.6p for A; ^ 1. The proof of Lemma l2.ll is 
concluded. □ 

Remark 2.2. Arguing as in [8], Lemma 3.2, it is possible to avoid the passage to the limit 
in the proof of Lemma 12.11 obtaining therefore a function h, solution to the homogeneous 
equation V • a{X, V/i) = 0, for the original vector field a. For our purposes though, it suffices 
to obtain an equation within the same universal class of a, p.6p . 

3 Optimal BMO estimates 

In this section we shall establish optimal a priori estimates for solutions to 

-V-a{X,Du) = f€Ll^^^{B^), 

which corresponds to the lower borderline integrability condition on /. In particular, 
bounds cannot be achieved under such a weak hypothesis. We recall that a measurable 
function / is said to belong to the weak-LP(i?i) space, denoted by L^^,^^{Bi), if there exists 
a constant K > for which 

L^{{XGB,:\f{X)\>T})<^. (3.1) 



The infimum of all > for which (|3.ip holds is defined to be the weak-L^ norm of u and it 
is denoted by (Bt)- Weak spaces play a fundamental role in Harmonic Analysis, 

weak ^ ^ ^ 

in particular in the theory of singular integrals. It is well known that ^ ^wcak' ^^^O' if 
A4 denotes the Hardy-Littlewood maximal operator, then M.{f) € -Z^^cak provided f £ L^, 
and such a result is optimal in the sense that Ai{f) may not belong to L^. This is the main 
reason for which Calderon-Zygmung theory fails for sources in L^. 

To motivate the result of this section, we invite the readers to notice that a careful 
inference in the kernel from ()1.3p revels a lower borderline condition for the source function 
f. In fact, r S L*^ for any r < but F L"-2 . That is, by Holder inequality, 

Nf G L°° whenever f £ L^'^", 

since ^ is the dual exponent of When / G , n > 3, Calderon-Zygmund estimate 

(|1.4p reveals that 

u e ly^'t ^ W^'"" ^ LI, (3.2) 
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for any 1 < g < oo. That is, it provides an almost optimal regularity result. By a duality 
argument, one finds out that it is impossible to bound the LJ^^-norm of u by the L2" norm 
of /. However, an application of Poincare inequality combined with (j3.2p gives 



\u - Ur^dX <Cn [ \Vu\''dX < C„||/||^n , (3.3) 

J Br 



where, Ur denotes the mean of u over Br, i.e., Ur := j-^ ^^^^^^ udY. 

Recall a function u £ L^{Bi) for which there exists a constant K > such that 

4 \u-Ur\<K, (3.4) 

JBriXo) 

for every Xq G Bi and < r < dist(Xo, (9-Bi), is said to belong to the BMO space. The 
infimum of all K > for which ()3.4p holds is defined to be the BMO-norm of u and it is 
denoted by ||n||BMO- 

The BMO space was originally introduced by John and Nirenberg in |13j. In that very 
same paper, John and Nirenberg proved the following fundamental estimate: if ||it||BMO ^ 
then there exist positive dimensional constants a and /3 such that 

/ e°l"-"il(iX < ^. (3.5) 
Jbi 

The original motivation for studying these functions apparently came from the theory of 
elasticity, [12]. Interestingly enough, John-Nirenberg's estimate for BMO functions (j3.5p 
is used by Moser as a key ingredient in his striking proof of Harnack inequality for diver- 
gence form uniform elliptic equations. Both Jonh-Nirenberg and Moser works were published 
simultaneously in the same issue: Comm. Pure Appl. Math. Vol XIV, back in 1961. 

Through the years, BMO space and its analogues have been shown to enjoy many other 
properties, with deep applications in analysis. For our purposes, it is elucidative to think the 
BMO space as the correct substitute for as the endpoint of the spaces as p f +00. 

In what follows, we will establish the corresponding sharp BMO estimate for solutions to 
p-degenerate elliptic equations 

-V-a(X,Dn) = /GL3^^,(Si). (3.6) 
where a satisfies the standard structural condition ()1.6p . 
Theorem 3.1. Let u G W^'^{Bi) he a solution to 

-V ■a{X,Du) = f{X). 



Assume a satisfies (II. 6p and f G ^^eafc(-^i)- Then u G BM0{Bi/2)- Furth 



ermore, 



1 



\\u\\bMO{B^/2) ^ C 11/11 "n' + 

V ^Lak(Bl) ^ 

for a constant C that depends only on n p, A and A. 
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In view of the parallel described above to the linear theory, the estimate from Theorem 
13.11 should be optimal. Indeed, this is the case. For instance, say, for p < n, if we set 

f{X) = \X\~P, it is easy to see that / G -^^eak- Solving ApU = f with constant boundary 
data on dBi one finds u{X) = Cn,p ■ In \X\, which is in BMO but not in L°°. 

The proof of Theorem 13.11 will be based on the compactness result granted in Lemma 12.11 
and an iterative scheme. Next Lemma is pivotal to our strategy. 

Lemma 3.2. Let u G W^'P{Bi) be a weak solution to (jl.Sp . with j- \u\PdX < 1. There 
exist constants 0<eo^l) 0<Ao^l/2, that depend only on n, p, A and X, such that if 



then 



4 \u{X) - (/ udY)\PdX < 1. 



Proof. Initially let us recall a general inequality: 



-I U--1 udY dX<2P-[ \u-j\PdX, 

J Br J Br J Br 



(3.7) 
(3.8) 

(3.9) 



for any u G and any real number 7. Indeed, by triangular inequality, 

i/p 



i - ifs 



\u — 7| 



udY — 7 



< 



< 2 



''dx) + / 

-I \u--ffdXj +-[ \u--f\dX 

J Bf / J Bf 

(-[ \u--f\PdX 



In view of Lemma |2. 11 let /i be a solution to the homogeneous equation in B1/2 such that 

/ \u{X)-h{XrdX<-^, (3.10) 



for Ao ^ 1/2 to be regulated soon. Such a choice will determine Eq. Notice that (I3.10p 

implies \h{X)\PdX < C, thus, by regularity theory for homogeneous equation, there 

1/2 

exists a constant C > universal such that 

\h{X) - h{0)\ < C7|Xpo, 
where C that depends only on n, p, A and A. Next, for Aq ^ 1/2 to be chosen, we estimate 

/ \u{X)-h{0)\PdX < 2P'U-[ \u{X) -h{X)\PdX +-[ \h{X)-hiO)\Pdx\ 



< 



9 • 2P 



+ C2P~' • a: 



1 \O!0P 



(3.11) 
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Now we can choose Aq, depending on dimension n and p, A and A so small that 

■ A"«^ < (3.12) 

and the proof of Lemma 13.21 follows from (j3.1ip and (j3.9p . □ 

Proof of Theorem \3.1[ Let u be a weak solution to 

-V ■a{X,Du) = f{X), inBi. 

The proof starts off with a renormalization. Let Eq be the universal constant from Lemma 
If we change u by nu, with k ^ 1, so small that 



and 



/ \K-u\PdX 

Jbi 



< 1, 



we can assume u and / are under the hypotheses of Lemma |3.2[ In the sequel, we will show 



\u-CkfdX<l, V/cGN. 



(3.13) 



^0 



Here Aq is the universal number from Lemma 13.21 and Ck denotes the average of u over the 
ball of radio Ag, i.e., 

:= / u{X)dX. 



Ck 



•^0 



We show (j3.13|) by induction. The case k = 1 follows directly from Lemma |3.2[ Assume we 
have verified p.l3p for k. We define the real function u : i^i — >■ M by 

v{X) := n(AgX) - Ck 



We also define 



a^k{X,C) := a{X'^X,C) and 



(3.14) 



(3.15) 



Notice that a^^k is also p-degenerate elliptic, with the same ellipticity constants as a. From 
the induction assumption, we have 



-I \v{X)\PdX = i |u(y) -Cfc|W < 1. 



Easily one verifies that 



If we label f^k := Aq'^ [/(AqA")], a direct computation reveals 

£"({XGSi:|/,g|>T}) = ({X G i?i : I/I > • Ao 



L^MX) 



< A: 



kp 



a.e. in B\. 



(3.16) 



(3.17) 



nk 



(3.18) 



< 



weak 



TV 
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That is 



fx^^Wj <\\f\\j <eo. (3.19) 

weak weak 



(3.20) 



We have verified that v is under the hypotheses of Lemma l3.2| which assures 
-[ \v{X) - d vdY)\PdX = -I \u{X) - Ck+i)\^dX < 1. 

This concludes the proof of (j3.13p . Finally, given < r ^ 1, let m G N be such that 

If we label Ur ■= j-^ udY, we estimate, 

■f \u-Ur\PdX < 2P-1 \u-XoUr\PdX 

J Bj, J Br 

< ^4 \u-Cm\''dX 

< c. 

The proof of Theorem 13. II can be now concluded my means of a standard covering argument, 
which we shall omit here. □ 

We finish up this section by highlighting once more that the strategy used in our reasoning 
to establish Theorem 13. II is indeed quite flexible. It is based on a fine scaling balance between 
the norm of the source / and the homogeneity of the equation itself. This indicates that 
similar analysis should be possible to be carried on for equations with measure data, provided 
the solution already lies in a proper Sobolev space, under the classical diffusion assumption 
\f\{Br) < Cr'^~P, for any ball B,. of radius r. For that, though, one needs to revisit the proof 
of Lemma 12.11 and work under appropriate notion of solutions through truncation. We do 
not intend to pursue that in this present paper. 

4 regularity 

In this section we turn our attention to optimal regularity estimates to Equation (jl.Sp when 

-+e 

the source function / lies in a slightly better space, say, / G -^^Xcak- ^^^^ case, heuris- 
tic scaling methods indicate that weak solutions should be locally bounded. Indeed, under 
slightly stronger assumptions on /, boundedness or even continuity of solutions can be deliv- 
ered by known methods, for instance through Serrin's Harnack inequality jl7j. Nevertheless 
this approach hardly reveals the sharp Holder exponent of continuity of the solution. 

In this section we still work under assumption (jl.6p . As we have already invoked, it is 
classical, see for instance |17] . that W^'^ solutions to the homogeneous equation 

- V • a{X, Du) = 0, in Bi, (4.1) 

are ao-Holder continuous in i?i/2 and 

lkllc"o{Bi/2) < C(n, A, A,p)||n||ip(Bj). (4.2) 
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The optimal exponent ao in (|4.2p depends only upon dimension, p and ellipticity constants 
A, and A. In general ao < 1 and its precise value is unknown. 

Theorem 4.1. Let u G W^''^{Bi) be a solution to 

-V ■a{X,Du) = f{X) (4.3) 

— 

Assume and f G LJ^^^{Bi), 1< 9 < p. Then u G C"(Si/2), for 

. r P - I _ 

a = mm{ • — )' (4-4) 

p — l a 

where ckq "is the universal optimal Holder exponent for solutions to —V • a{X, Du) = 0. 
Furthermore, 

IKIIc-iBi/a) ^ C{n,\K,p,e) { 11/11 + ) . 

The sharp relation in (j4.4p should be read as follows: 

If ^•Y<«o then nGCig^'"^. (45) 
If • > ao then u G Cj^j., for any a < uq. 

The proof of Theorem 14.11 will be given in subsection 14.11 below. Optimality of the thesis 
of Theorem 14.11 can be checked directly by computing in the unit ball, Bi 

P e-i p e-^ 

A,\X\ — -—=c\X\-o^LJ^^. 

It is interesting to notice that l-'^l"* is not in the classical Lebesgue space L ' p . A valuable 
feature of Theorem 14.11 is the fact that it provides universal bounds, i.e.. Holder estimates 
that depend only on ellipticity and p-degeneracy feature of the operator. This is particularly 
important in homogenization problems. However, under continuity (or some sort of VMO 
condition) on the medium, we can show that solutions to the homonegenous equation 

-V • a{X,Du) = 0, 

are C° for every a < 1. Indeed this fact is an immediate consequence of our next Theorem. 

In the sequel, we shall slightly improve the thesis of Theorem 14.11 provided the medium 
has some sort of continuity property. For simplicity purposes, for the next Theorem, we shall 
work under classical continuity assumption on the operator a with respect to the X variable. 
That is, there exists a modulus of continuity r such that 

p-i 



\aiX,0-aiY,()\<T{\X-Y\Mr'- (C) 



We remark that under the structural assumption (jl.6p solutions to the homogeneous, 
constant coefficient equation have a priori C^'*^ estimates for Xq G -B1/2 fixed. That is 

-Va{Xo,Dh) = 0, Bi implies \\h\\ci,.^By,) < C{n,p, X, A)\\h\\LPiB,), (4.6) 

for some < e < 1 that depends only on p, n, A and A, see, for instance, [7]. 



11 



Theorem 4.2. Let u G W^'^{Bi) be a solution to 



-V-a{X,Du) = f{X). (4.7) 

n n p 6 — 1 

Assume ([L6]), © and that f e L p{Bi), 1 < 9 < p. Then u £ Cp-^' « (^1/2) and 
furthermore, 

\M p.e-i^ < C(n, A, A,p,T, (9) I 11/11 e.^ + II'^IIlp(Si) ) • 



Before delivering the proofs of Theorem 14.11 and Theorem 14.21 let us make few comments 
about our regularity estimates. Initially, as in Theorem 13. H it seems reasonable to 

On — 'p 

establish the same optimal result for measure data /, provided \ f\{Br) < Cr » , for any 
ball of radius r. As for Theorem 14. 2t continuity condition can be greatly relaxed. In fact 
all we need is a sort of Cordes-Nirenberg type of condition: there exists a universal constant 
(5* > such that 

ia(x,e)-a(o,e)i <5*ier-^ 

The upper threshold case for continuity theory, / G L", is a delicate issue, see [H]. At this 
point, though, an interesting consequence of Theorem 14.11 is that solutions to 

-Va(X,Z)u) = /GL:,,k(^i), 

for measurable coefficients equations, has almost the same modulus of continuity as a- 
harmonic functions, i.e., solutions to —Va{X^Dh) = 0. That is, if a-harmonic functions 
in Bi are locally C"o, then solutions to -Va{X,Du) = / G -^^wcak(-^i) are locally , for 
any < /3 < oq- The same analysis employed in Theorem 14.21 gives that for equations with 
continuous coefficients, solutions to —Va{X,Du) = / S -^weak(^i) locally C^, for any 
< /? < 1. 



4.1 Proof of Theorem 14.11 

We revisit the proof of Lemma 13.21 Suppose ;f^^ < 1 and for q = 9 ■ p/n, 

Sl>\\f\\L^ JBO>Cn||/|U , 

weak \ ^ / J^P 

weak 

with ei > to be chosen. Prom Lemma |2 . 1 1 there exists a function h, solution to 

-V • a(X, D/i) = 0, inBi/2 

such that 

/ \u{X) - h{X)\P < 5i. 

The latter choice for 5i determines ei through the compactness Lemma |2. 11 Since ||/i||Lp is 
under control, the regularity theory for homogeneous equation assures h G C°"(i?i/3) and 
for a universal constant C > 0, 

\h{X) - /i(0)| < C\XY'\ 
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We can readily estimate 

\u{X) - h{0)\PdX < 2^"^ \u{X) -h{X)\PdX +-[ \h{X) -h{0)\PdX 



-1 \PO() 



(4.8) 



Now, fixed a < ao we can choose Ai ^ 1 universally small so that 



1 - 1 y J 

Once Ai is chosen as indicated above, we select 61 (and therefore ei) as 

= ^A'^+". (4.10) 
If we combine (|4.8|) . (|4.9|) and (|4.1U|) we conclude that 

/ \u{X) - h{0)\PdX < Af , (4.11) 
Jb^^ 

provided 



weak ^ ' 

for < ei ^ 1 that depends only on dimension, p X, A and a < oq. In addition, from the 
regularity theory for homogeneous equation, 

IMO)|<C, (4.13) 

for a universal constant C > 0. 

We remind that the assumptions -f „ < 1 and (b.) < can be reached 

-^1 weak ^ ^ ^ 

by a simple change of scaling and normalization. Thus, with no loss of generality, we can 
work under these hypotheses. 

In the sequel we shall prove that there exists a convergent sequence {^fcjfcgN C M for 
which 

|n(X) -/ifclW < A^^". (4.14) 



^1 



As before, we will verify (I4.14p by induction. The case A; = 1 is precisely (14. lip , with 
Hi = h{0). Suppose we have checked (|4.14|) for /c = 1, 2, • • • , m. Define 

»m := "'""g^-"" (4.15) 

With the same notation as in (I3.15p . we readily verify, as in (I3.17p . that 

\Lx^MX)\ < A^'^-^^-')"! l/(Ar^)| =: fm{X). (4.16) 
One easily estimates, for any r > 0, 

£«({XGi?i:|/^|>r}) = f £ i?i : I/I > ) " 



-m-n 



^ weak 



yiii\^li-yp-Llui\q ^ 



(4.17) 
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in view of the sharp assumption (j4.4p . We have shown that v is entitled to the conclusion in 
(|4.1ip . Let hm be the solution to the homogeneous problem that is -^/^-close to v in -B1/2 
in the L^-distance. We label hra{0) = tm and, as in ()4.13p . \ tm\ < C for a universal constant. 
Applying (|4.1ip to v we find 

/ \v{X) -tml^'dX < XP". (4.18) 

Rescaling (|4.18p back yields 

HX) - i^^r,^ + XrtrnWdX < Af ("^+^). (4.19) 



i 

^1 



Therefore, the induction step for (|4.14p is verified by taking 
Indeed {nk}keN is a convergent sequence, because we estimate 



as A; — 7- oo. Finally, if we define 



fi := lim 



and < r < 1 is arbitrary, estimate (14.14P gives 

-j \uiX) - -p\PdX < CrP'^; 

therefore u is a-Holder continuous at the origin. The proof of Theorem 14.11 follows now via 
standard covering arguments, which we omit here. □ 

4.2 Proof of Theorem [4721 

For convenience, let us label q := 6 ■ ^ > ^. The proof of Theorem 14.21 is based on the 
following refinement of the Compactness Lemma |2.1[ 

Lemma 4.3. Let u € VF^'^(i?i) be a weak solution to (jl.Sp . with f \u\PdX < 1. Given 
5 > 0, there exists a < e ^ 1, depending on depending only on p, n, X, A and 6 such that 
if 

Wfh^ m^)<^, <^^d \a{X,i)-a{Q,i)\<e\i\P-\ (4.20) 
then there exists a function h in B1/2 solution to 

- V • a{Dh) = 0, in B1/2, (4.21) 

for some constant coefficient vector field a satisfying ()1.6p with the same ellipticity constants 
X and A, such that 

\u{X) - h{X)\PdX < 6P. 



i 

•^^1/2 
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Proof. As before, let us assume, searching for a contradiction, that the thesis of the Lemma 
fails. If so, there would exist a 6q > and sequences 

UkGW'^^iBi), and A G L^^^^l^i), 

with 

\uk{X)\'PdX < 1, (4.22) 



JBj 

for all A; > 1, 

-V-ak{X,Duk) = fkinBi, (4.23) 

where Ok satisfies ()1.6p . and 

IIAIUnSi) + \^\'-''\ak{X,0 - afc(0,e)| = o(l), (4.24) 

as A; — )■ 0; however 

/ \ukiX)-h{X)\PdX>5o, (4.25) 

for any solution /i to a homogeneous, constant coefficient equation (j4.2ip . in i?i/2 ai^d all 
k > 1. Reasoning as indicated in the proof of Lemma |2.H we have 

\Vuk\^dX < C, 

for all A; > 1. Thus, up to a subsequence, there exists a function u G ^^^'^(^1/2) for which 

Uk ^ u in W^'^{Bi/2), Uk —f- u in L^{Bi/2), and Vnfc(X) — \7u{X) a.e. in -Bi/2- 

(4.26) 

Also, by Ascoli Theorem, there exists a subsequence under which 0^^.(0, •) — )• a(0, •) locally 
uniformly. Thus, for any X G -B1/2, 

|afc^,(X,e)- 0(0,01 < |afc^.(X,e)-afc,(0,OI + |afc,(0,0-a(0,OI =0(1), (4.27) 

that is, akj{X,-) — )• a(0, ■) locally uniformly. Finally, given a test function cf) G Wq'^{Bi/2), 
in view of (0231), (t^^ and (lOTll we have 

/ ak{X,Vuk)-V(pdX = I fkifdX 

/ a(0,Vn) • V(/)(iX + o(l), 

as A: —7- 00. Since (j) was arbitrary, we conclude u is a solution to a constant coefficient 
equation in -Bi/2- Finally we reach a contradiction in (j4.25p for A; ^ 1. □ 



The main difference between Lemma 12.11 and Lemma 14.31 is the fact that the former 
provides existence of a C"^ function close to u under smallness assumptions on the data. The 
latter gives a function near u under smallness assumptions that also involve continuity of 
the medium. Thus, the following version of Lemma 13.21 can be proven by similar arguments 
used to establish estimate ()4.1ip . 
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Lemma 4.4. Let u G W^'^{Bi) be a weak solution to (|1.5p . with f \uY'dX < 1. Given 
a < 1, there exist constants < Eq <^ 1, < Aq ^ 1/2 and cq G M such that if 



< ^0 and |a(X,e)-a(0,e)| <eo|er\ (4-28) 

then 

\u{X) - co\PdX < Xl"". (4.29) 



Proof. For 5 > to be regulated a posteriori, let h he a solution to a constant coefficient 
equation assured by Lemma |4.3[ that is 5-close to u in the L^-norm. From C^'*^ regularity 
theory for constant coefficient equations, (|4.6p . there exists a constant C depending only on 
n,p, A and A such that 

\h{X) - h{0)\ < C\X\. 

Since ||/i||LP < C, by L°° bounds, 

\hm < a 

We now estimate, 

/ \u{X) - h{0)\PdX < 2P-m/ \h{X) - h{0)\PdX + -[ \u{X) - h{X)\Pdx] 

Since < a < 1, it is possible to select Aq small enough as to assure 

2^-^CAg < Ia^^ 

Once selected Aq, we set 

1 

which determines the smallness condition eq through the compactness Lemma |4.3[ □ 

Finally, the proof of Theorem 14.21 follows by the induction argument from Section 14. H 
having Lemma 14.41 as its starting basis. We omit the details here. 
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